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Shear thickening is a widespread phenomenon in suspension flow that, despite sustained study, is 
still the subject of much debate. The longstanding view that shear thickening is due to hydrody¬ 
namic clusters has been challenged by recent theory and simulations suggesting that contact forces 
dominate, not only in discontinuous, but also in continuous shear thickening. Here, we settle this 
dispute using shear reversal experiments on micron-sized silica and latex colloidal particles to mea¬ 
sure directly the hydrodynamic and contact force contributions to shear thickening. We find that 
contact forces dominate even continuous shear thickening. Computer simulations show that these 
forces most likely arise from frictional interactions. 


Shear thickening, the increase of viscosity with shear 
rate, is ubiquitous in complex fluids urn- In particular, 
it plays a controlling role in the flow of concentrated sus¬ 
pensions of micron-sized particles [5]. Such suspensions 
occur widely in applications, from ceramics and bullet¬ 
proof armor, through cement and even chocolate. Shear 
thickening has also become a mainstay of popular science, 
in the form of running on a pool of corn starch solution. 

Despite sustained study, the mechanism of suspension 
shear thickening is still disputed. A long-standing view is 
that thickening, especially if it is continuous with shear 
rate, is predominantly driven by hydrodynamic interac¬ 
tions laniisiiT]. In some theories, the increase in viscos¬ 
ity arises directly from enhanced dissipation in the nar¬ 
row lubrication films between particles mil]; in others, 
it is caused by the formation of transient particle clus¬ 
ters (“hydroclusters”) msiiaiaiini. In support, exper¬ 
iments have identified putative hydroclusters mmm 
and found an increase in hydrodynamic stresses during 
thickening [IMS]. 

Recently, contact forces have been shown to mediate a 
discontinuous jump in viscosity with shear rate in dense 
suspensions of non-Brownian particles at volume frac¬ 
tions greater than random loose packing 0 > 0.58 HHH- 
uni- More controversially, simulations Enm and the¬ 
ories [23j (24] suggest that direct frictional contact can 
also lead to continuous shear thickening in moderately 
dense colloids with (f < 0.58. This proposal has not yet 
gained wide acceptance: conceptually, it seems harder 
to find a role for such forces without the formation of 
system-spanning contact networks [4]. 

Experiments that can dissect the relative contributions 
of hydrodynamics and contact stresses can settle this is¬ 
sue definitively. Here, we demonstrate that shear reversal 
techniques, pioneered by Gadala-Maria et al. [25j|26], can 
distinguish the relative contributions of these stresses in 
non-inertial shear thickening systems. The basic idea is 
simple but powerful: immediately upon shear reversal, 
instantaneous contact stresses vanish, but hydrodynamic 
stresses remain because of Stokes flow reversibility. Thus, 


monitoring time-dependent stresses after reversal will re¬ 
veal the relative roles of these two interactions in thick¬ 
ening. 

We do so in two canonical colloidal systems: charge- 
stabilized silica (Seikisui Ghemical) with diameter 2.0 fim 
suspended in a mixture of glycerol and water (viscosity 
rjo = 0.98 Pa.s at 20°G) at volume fraction (f = 0.49, 
and polymethylmethacrylate (PMMA) particles with di¬ 
ameter 1.4 fim sterically stabilized with poly-dimethyl- 
diphenyl siloxane (PDV-2335, Gelest) with chain length 

50 nm m suspended in PDV-2331 (Gelest, viscosity 
rjo = 1.78 Pa.s at 20 °G ) at 0 = 0.51. 

Shear reversal measurements were performed in an 
ARES strain-controlled rheometer (Rheometric Scien¬ 
tific) with roughened cone-plate geometry (25 mm, 0.1 
rad) modified with a DAQ directly connected to the ana¬ 
logue output of the stress and strain sensors [28] . Directly 
measuring the output of these sensors bypasses signal 
processing that hinders the instantaneous measurement 
of the system’s short-time response. 

In our protocol, figure [^a), a positive shear rate 7 
is imposed until the accumulated strain 7 reaches 10 , 
after which the shear is reversed and a negative shear 
rate is imposed to accumulate the same amount of strain. 
Just before shear reversal, the suspension is in a steady 
state. Immediately upon reversal, the structure is un¬ 
changed, so that hydrodynamic forces will remain iden¬ 
tical in magnitude but reversed in direction. In con¬ 
trast, contact forces in hard-sphere systems will drop 
to zero immediately after reversal, figure Bb). This 
qualitative difference between these two forces allows us 
to disentangle their separate contributions by measur¬ 
ing the transient stress upon reversal. Eigure l^a) illus¬ 
trates schematically the instantaneous relative viscosity 
r]{t) = cr(t)/(7(t)7o) as a function of time (or, equiva¬ 
lently, accumulated strain) of a purely Newtonian fluid 
(green line) and the response where contact forces con¬ 
tribute to the total stress (red line). (See Supplemental 
Information for calibration.) 

A rate-sweep measurement of the charge-stabilized sil- 
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FIG. 1 . (a) Schematic of our shear reversal protocol. The ap¬ 

plied strain 7 is a triangle wave with peak-to-peak amplitude 
of 10 . The applied strain rate 7(t) is a square wave whose 
amplitude we vary in our experiments. The instantaneous 
stress response, cr, of a viscous fluid and a contact-dominated 
fluid are shown in the green and red, respectively. (The same 
color scheme applies throughout all parts of this hgure.) The 
instantaneous relative viscosity r]{t) = cr(t)/(7(t)?7o). (b) 

Schematic illustrating the difference between hydrodynamic 
forces, which stay the same in magnitude immediately upon 
reversal, and contact forces, which drop to zero immediately 
upon reversal. 


ica suspension gives its flow curve, 77(7), figure[^a) (thick 
grey line), which shows shear thickening at 7 > 0.3 s“^. 
Next, we monitor the viscosity after shear reversal, r^rev, 
as a function of accumulated strain after reversal, 7, at 
four representative shear rates 7 = 0.50,1.00,2.00 and 
5.00 s“^, in the shear thickening regime, figure [^b). Im¬ 
mediately before reversal, 77 has a value (not shown) equal 
to that measured in the steady state (7 = 10). The vis¬ 
cosity immediately after reversal drops to a value, 77 !?^^, 
that remains constant, to within experimental error, until 
7 > 0.1. We take 77^^^ to be the hydrodynamic contri¬ 
bution to the total steady-state viscosity. From 7 ^ 0.3, 
77rev rises, reaching a steady-state value, 77 ^^, that is the 
same as the steady-state viscosity before reversal. We 
take 77 ^^ — 77 j?g^ to be the contact contribution. 

The hydrodynamic and contact contributions to the 
total viscosity so obtained are plotted as a function of 7 
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FIG. 2. Relative viscosity, 77, vs strain rate, 7, for the (a) 
silica and (b) PMMA suspensions. Shown are rate sweep 
measurements (thick grey line), steady state relative viscosi¬ 
ties from flow reversal (red) decomposed into contributions 
from hydrodynamic (blue) and contact (green) interactions, 
(b) and (d): Instantaneous relative viscosity after shear rever¬ 
sal, Tjrevi'y) = cr( 7 )/( 777 o), VS Strain, 7, at different applied 7. 
The hydrodynamic component of the relative viscosity, r]rev^ 
is taken to be the average of the relative viscosity over the 
range 0.01 < 7 < 0 . 2 . The contact component is taken to be 
the difference between the steady state value at 7 = 10, 77^^, 
and 77 °e^;. 


in figure [^a) . Strikingly, while the contact contribution 
increases with 7 , the hydrodynamic contribution remains 
constant as the suspension shear thickens. This demon¬ 
strates the essential role played by contact forces in the 
continuous shear thickening of this silica suspension. 

Contact forces depend sensitively on the nature of par¬ 
ticle surfaces. We therefore repeated our experiments 
using a PMMA suspension. In contrast to our silica 
particles, which are charge stabilized, these PMMA par¬ 
ticles are stabilized by (neutral) surface ‘hairs’. Our 
(j) = 0.51 PMMA suspension shear thickens more readily, 
at 7 0.01 s“^, figure ^c). Again using shear-reversal 

to measure figure ^d), and analyzing the data in 
the same way to obtain the hydrodynamic and contact 
contributions, figure [^c), we find essentially the same 
pattern of behavior already uncovered for the 0 = 0.49 
silica suspension. Thus, contact forces also dominate the 
continuous shear thickening in this sterically-stabilized 
system. 

To validate this interpretation, and to probe the nature 
of these contact forces, we conducted computer simula¬ 
tions of repulsive spheres (radius a) in which the short- 
range lubrication and repulsive contact forces were mim- 
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FIG. 3 . Critical load model simulations of particles with hydrodynamic and frictional interactions, (a) Relative viscosity 77(7) 
versus strain 7 at four different shear rates (given in the legend, see text for the definition of 70); compare experimental results 
in figure l^b) and (d). (b) Directly-calculated hydrodynamic and contact contributions to the total viscosity as functions of 
the accumulated strain after reversal for the case of 7 = O.7470. (c) The total viscosity as a function of shear rate, and its 
decomposition into hydrodynamic and contact contributions done in two ways: directly calculated from simulations, and ‘read 
off’ the data shown in part (a) in the same way as the experimental data is analyzed into these contributions, for which see 
figure l^b). 


icked using a previously-established critical load model 
[ 21 ], implemented in a classical discrete element method 
code (details of which are identical to those in |29|; also, 
see Supplemental Information). In this model, frictional 
interactions appear beyond a critical normal force be¬ 
tween particles, which also sets a shear rate scale, 

7o = /{67rr]oa‘^). We explored shear reversal at 

0 = 0.51. 

The observed evolution of the viscosity with accumu¬ 
lated strain after reversal, figure |^a), is qualitatively 
similar to experiments. In the simulations, however, we 
can follow the hydrodynamic and contact contributions 
as a function of accumulated reversed strain by direct 
evaluation; figure [^b) shows the results for 7 = 0 . 7400 . 
Consistent with the interpretation we have offered for 
our experimental data, the contact contribution drops 
essentially to zero immediately after reversal, and only 
increases back to its steady-state magnitude after an ac¬ 
cumulated stain of order unity. By contrast, the hydro- 
dynamic contribution remains more or less constant. 

We cannot access directly the hydrodynamic and con¬ 
tact stresses in our experiments. Instead, we have to 
‘read off’ these contributions to the reversed viscosity 
from the data as 77 ^^^ and 77 ^^ — 77 ^^^, respectively, fig¬ 
ure Hb). We validate this procedure using simulations. 
First, we ‘read off’ the hydrodynamic and contact con¬ 
tributions to the reversed viscosity data, figure [^a), in 
exactly the same way as in experiments, giving results, 
figure 3 c), analogous to the experimental data shown in 
figure 2 ^a) and (c). Secondly, we evaluated these contri¬ 
butions directly from the raw simulation data, and over¬ 
laid the results in the same graph, figure |^c). The near 
identity of the two sets of results validates our identi¬ 
fication of 77 ^g^ and 77 ^^ — 77 !?^^ with the hydrodynamic 
and contact contributions, respectively. Thus, simula¬ 
tions confirm the key role played by contact forces in 


continuous shear thickening. 

Traditionally, shear thickening is supposed to occur be¬ 
yond a critical onset strain rate | 30 ]. Recent theories, 
simulations and experiments 0 El Ell Ea El suggest 
that, instead, the phenomenon appears at a critical on¬ 
set shear stress. Indeed, this assumption has been built 
into the critical load model used in our simulations. To 
explore the role of stress in the onset of shear thickening, 
we performed rate sweep measurements in our two sys¬ 
tems at different temperatures, T, to change the solvent 
viscosity, 770 (F), over an order of magnitude. 

We find that as the solvent viscosity varies, the onset of 
thickening (dashed lines in figure |^a) and (c)) occurs at 
different strain rates for both of our suspensions. Thus, 
in the silica suspension, figure |^a) , the highest solvent 
viscosity (lowest temperature) data set shows shear thick¬ 
ening at 7 > 200 s“^, but this onset progressively moves 
to higher rates as we lower the solvent viscosity, until 
in the lowest viscosity data set, we barely see thickening 
over our range of shear rates, but observe a small de¬ 
gree of shear thinning instead. Similarly, in the PMMA 
suspension, figure [^c) , we see the onset of thickening at 
around 0.05 s~^ in the highest viscosity data set; this 
onset also moves to higher shear rates as the solvent vis¬ 
cosity decreases. 

Each of these two data sets, however, can be scaled 
onto a master curve if we plot the relative viscosity 
77(7) = (7/(7770 (T)) against the shear stress a [H]. The 
master curve for the silica particles, figure [^b), indeed 
shows a single onset stress at ~ 10 Pa, while that for 
the PMMA system, figure |^d) , shows an onset stress at 
1 Pa. Thus, in both systems, our experiments show 
that shear thickening occurs above a critical stress. The 
difference in the magnitude of the onset stress of the two 
suspensions is consistent with the expectation that the 
load needed to press particles into contact should be sen- 
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sitive to details of the stabilization mechanism [32]. 

Our observation of contact-driven shear thickening 
contradicts the conclusions drawn from rheo-SANS [T4| 
and scattering dichroism measurements on Brownian 
hard spheres m- The rheo-SANS data provides an ap¬ 
proximation of the hydrodynamic stresses through a non¬ 
trivial analysis of the structural anisotropy in the suspen¬ 
sion. However, even the latest rheo-SANS measurements 
show thickening that is only a small fraction of that ob¬ 
served in rheometry m- In the scattering dichroism 
measurements, the total suspension viscosity, 77 , and the 
viscosity contribution from Brownian stresses, rjB^ can 
be measured independently. The authors found that the 
contribution to the viscosity from other stresses, 77 — 77 ^, 
increased upon shear thickening, and attributed this to 
an increase in hydrodynamic stresses m- However, the 
remaining viscosity may in fact be dominated by contact 
contributions, an interpretation that is in alignment with 
results from dense non-Brownian suspensions [33] as well 
as the data presented in the current work. 

Scattering dichroism measurements [34] also found ev¬ 
idence of particle clusters in the shear-thickened state 
(later confirmed by real-space imaging [T]), which were 
supposed to originate from large lubrication stresses be¬ 
tween particles. The origin of these clusters remains un¬ 
clear, however - the present work suggests that any clus¬ 
tering is likely to involve frictional contacts. 

To conclude, we have used shear reversal to show that, 
for two canonical model colloidal systems, continuous 
shear thickening does not originate from hydrodynamic 
interactions but from the formation of particle contacts. 
The onset of contact formation, and thus shear thick¬ 
ening, is found to occur above a critical stress whose 
value is sensitive to whether the particles are charge- or 
sterically-stabilized. Further work is needed to estab¬ 
lish whether shear thickening of nanoparticles, for which 
Brownian stresses are not negligible, is also driven by 
contact formation; this should be possible with shear re¬ 
versal. 

While shear reversal alone provides no a priori infor¬ 
mation about the nature of the interactions between con¬ 
tacting particles, the quantitative agreement between the 
experimental and simulation viscosity response after re¬ 
versal, figure 2(b), (d) and 3(a), strongly suggests that 
static friction is present. Indeed, evidence for frictional 
contacts was recently found in a similar model system [5] . 
A careful study of the effect of particle surface properties 
- specifically surface roughness - on suspension rheology 
would be required to validate this claim, which in turn 
motivates the need for robust measurements of particle 
friction. 
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FIG. 4. Scaling of the relative viscosity with shear stress, 
(a) and (c): Viscosity versus shear rate for silica (0 = 0.43) 
and PMMA (0 = 0.46) suspensions at different temperatures. 
Solvent viscosities and temperatures are as labeled in (b) and 
M) (for silica, we used a different solvent composition to figure 
[H see Supplementary Information). Due to the temperature 
dependence of the solvent viscosity the initial suspension vis¬ 
cosity and onset strain rate (arrows) of the thickening are 
substantially altered. (b,d) Relative viscosity vs stress show¬ 
ing collapse of the thickening data onto a master curve. 
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